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Abstract 

In disordered vortices in p-wave superconductors the two new random-matrix ensembles 
may be realized: B and Dill-odd (of so(2N + 1) and so(4JV + 2)/u(27V + 1) matrices respec- 
tively). We predict these ensembles from an explicit analysis of the symmetries of Bogoliubov- 
deGennes equations in three examples of vortices with different p-wave order parameters. A 
characteristic feature of the novel symmetry classes is a quasiparticle level at zero energy. Class 
B is realized when the time-reversal symmetry is broken, and class Dill-odd when the time- 
reversal symmetry is preserved. We also suggest that the main contribution to disordering the 
vortex spectrum comes from the distortion of the order parameter around impurities. 

Since Wigner's modeling Hamiltonians of complex nuclei by random matrices Q, the random- 
matrix theory (RMT) has played an important role in studying mesoscopic systems. In many cases, 
a chaotic (non-integrable) mesoscopic system may be accurately described by RMT. The super- 
symmetric technique by Efetov provides a microscopic explanation of the RMT approximation 
for disordered systems. 

According to the RMT approximation, the only characteristics of the system affecting the 
eigenvalue correlations at small energy scales are its symmetries. Therefore, a problem arises of 
classifying symmetries of random-matrix ensembles. It has been suggested by several authors || 
[]]. ^ that random-matrix ensembles may be classified as corresponding symmetric spaces. The 
symmetric spaces may be divided into twelve infinite series reviewed in Table 1 (we split class .Dili 
in two subclasses: Dill-even and -Dill-odd) 0. Each of the symmetry classes may occur in one of 
the three forms: positive-curvature, negative-curvature and flat The corresponding Jacobians 
in the matrix space are expressed in terms of trigonometric, hyperbolic, and polynomial functions 
respectively. In the present paper we shall only discuss the RMT for Hamiltonians forming a linear 
space and therefore described by zero-curvature (flat) versions of RMT. 

The simplest examples of the RMT symmetry classes are the three Wigner-Dyson classes: 
unitary, orthogonal and symplectic (A, Al, and All, respectively, in Cartan's notation) ||. In 
these classes, the energy level correlations are invariant under translations in energies (at energy 
scales much smaller than the spectrum width), and the joint probability distribution of the energy 
levels u>i is 
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The parameter (3 determines the strength of level repulsion and takes values 2, 1, and 4 for the 
unitary, orthogonal, and symplectic ensembles, respectively. 
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Tabic 1: Symmetric spaces and universality classes of random- matrix ensembles: (a) the three 
Wigner-Dyson classes A (unitary), Al (orthogonal), All (symplectic) ||; (b) chiral classes Alll 
(unitary), BDI (orthogonal), CII (symplectic), with dimensions p < q have q — p zero modes ||; 
(c) superconducting classes C, D, CI, fill-even ||; (d) p-w&ve vortex classes Dill-odd, B have 
one zero mode 
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Other symmetry classes appear when there exists an additional symmetry relating energies E 
and —E. In the corresponding random-matrix ensembles, the levels LUi repel not only each other, 
but also their mirror images —u>i, and the Jacobian has the form 

It is now characterized by the two parameters f3 and a (the latter responsible for suppressing the 
density of states near E = 0). 

The three chiral classes (AIII, BDI, and CII) appear in systems with a chiral symmetry anti- 
commuting with the Hamiltonian 0. Four more classes (C, D, CI, and Dili) describe mesoscopic 
superconducting systems (with or without spin-rotational and time- reversal symmetries) Q . The 
present paper is devoted to the last two lines in Table 1: classes Dill-odd and B. We shall 
demonstrate that these two classes correspond to the symmetries of a p-w&ve vortex with or with- 
out time-reversal symmetry, respectively (the symmetry class B in p-wave vortices has also been 
predicted in |[o| ; a partial account of the present work appeared as a preprint jll| ) . 

Below we consider three particular p-wave vortices and describe their symmetries: the single- 
quantum vortex in the A-phase with a fixed spin orientation of the order parameter, the half- 
quantum vortex in the A-phase with rotating orientation of the order parameter, and the spin 
(fiuxless) vortex in the -B-phase. We explicitly demonstrate that in the first two examples the 
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symmetries of the Hamiltonian are of type B, and in the last example - of type 15111-odd. Then 
we briefly discuss microscopic requirements for the RMT limit and show that the main contribution 
to the level mixing comes from inhomogeneous suppression of the order parameter by impurities. 
Our estimates suggest that the RMT limit may be achieved in a moderately clean superconductor. 



1 Single-Quantum Vortex 



In this section we consider a vortex in the order parameter similar to the A phase of 3 He. Namely, 
we assume that the condensate wave function has the form 



i, ( ITT) + III) ) + «2y(|TT) - III) ) + d z ( |U) + III) 



(k x ± iky) , (3) 



and that the vector d defining the orientation of the triplet is fixed (without loss of generality, 
take it to be in the z direction. We also assume a fixed chirality of the order parameter and take 
positive sign in place of ± in (^) . Then the only allowed vortices are those involving a rotation of the 
phase of the order parameter, and the magnetic flux in such vortices is quantized in conventional 
superconducting flux quanta $0 = hc/2e. 

Below we analyze the low-lying energy levels obtained as solutions to Bogoliubov-deGennes 
equations for the mean-field Hamiltonian 



*-E»![ fe i^ + VM-.,]». + .t(A..& + A..£) 



A I a , a ^^j+h.c. , (4) 



where ^ a are the electron operators (a is the spin index), V(r) is the external potential of impuri- 
ties, A(r) is the electromagnetic vector potential, A x (r) and A y (v) are the coordinate-dependent 
components of the superconducting gap. [In the bulk, the preferred superconducting order is one 
of the two chiral components r\± — A x ± iA y , but in inhomogeneous systems, such as a vortex 
core, an admixture of the opposite component may be self-consistently generated $L2\ . We ac- 
count for this effect by allowing the two independent order parameters A x and A y .] Star (*) 
denotes the symmetrized ordering of the gradients and the order parameters A M [definition: 
A* B = (AB + BA)/2]. At infinity, the order parameter imposes the vortex boundary conditions: 

A x (r^oo,0) = A o e ± ^, A y (r -f oo, <j>) = , (5) 

where r and <f> are polar coordinates. Plus or minus signs in the exponent correspond to a positive 
or a negative single-quantum vortex. For an axially-symmetric vortex with the chirality of the order 
parameter non-self-consistently fixed (A y = iA x ), without the vector-potential A(r) and without 
disorder V(r), the low- lying eigenstates of the Hamiltonian (f|) have been found by Kopnin and 
Salomaa [Q. The spectrum is 

E n = nujQ , (p — wave) , n = 0, ±1, ±2, ... (6) 

with luq r~j A 2 /ef- This result should be compared with the spectrum of the vortex core in a 
s-wave superconductor [fl4|| : 

E n =(n + ^jujQ, (s-wave), n = 0, ±1, ±2, . . . (7) 
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The common feature of the spectra in the s-wave and p-wave cases is the symmetry with 
respect to zero energy. If we interpret holes in the negative-energy levels as excitations with 
positive energies (and with opposite spin), then this symmetry implies that the excitations are 
doubly degenerate in spin: to each spin-up excitation there corresponds a spin-down excitation 
at the same energy. For the s-wave vortex, this degeneracy is due to the full spin-rotation SU(2) 
symmetry. The p-wave Hamiltonian (^) has a reduced spin symmetry. Namely, it has the symmetry 
group 0(2) generated by rotations about the z-axis ("fj i— > e ,a $|, "Jj i— > e~ la ^> i) and by the spin 
flip v&j i ► ^> i, ^> i i ► vffj. This non-abelian group causes the two-fold degeneracy of all levels, 
except for the zero-energy level(s) where the symmetry 0(2) may mix the creation and annihilation 
operators for the same state. This symmetry is crucial for our discussion. Note that we have not 
included in the Hamiltonian neither the spin-orbit term (Ugo ■ [cr x p]), nor the Zeeman splitting 
H(r) ■ a. Either of these terms would break the spin symmetry 0(2), which would eventually result 
in a different universality class of the disordered system (type D with non-degenerate levels), if 
these terms are sufficiently strong. 

The difference between the s- and p-wave vortices is the zero-energy level in the p-wave case. 
It has been shown by Volovik that this level has a topological nature pjj| . Indeed, suppose we 
gradually increase disorder in the Hamiltonian (^). The levels shift and mix, but the degeneracy 
of the levels remains the same as long as the symmetry 0(2) is preserved. The total number of 
levels remains odd, and therefore the zero-energy level cannot shift if the final Hamiltonian is a 
continuous deformation of the original one (without disorder), i.e. if the topological class of the 
boundary conditions (^|) remains the same. 

Now we can identify the symmetry class of the Bogoliubov-deGennes Hamiltonian. The only 
symmetry of the mean-field Hamiltonian (|J) is the spin symmetry 0(2). The time-reversal symme- 
try is already broken by the vortex and by the pairing, and therefore neither the vector potential 
A(x) nor local deformations of A M can reduce the symmetry of the Hamiltonian. When projected 
onto spin-up excitations 7^ = J[u(r)4 r |(r) + v(r) 1 $> i(r)]d 2 r, the Bogoliubov-deGennes Hamiltonian 
for the two-component vector (u, v) takes the form: 




#BdG = I [a* , .„ , A* 



V(r) - e F 

HV X ) + g * (-iVy) 



& * Hv x ) + 1* * hv„ 

( " v+eA)2 + V(r) - e F 



(8) 



2m 

In an arbitrary orthonormal basis of electronic states, this Hamiltonian may be written as a matrix 

#BdG = y J • ( 9 ) 

From the hermiticity of the Hamiltonian, it follows that tv — h. From the explicit form of the 
p-wave pairing, A = — A T (it is here that the p-wave structure of the pairing is important; for 
s-wave pairing we would have A = A T instead) . These are the only restrictions on the Hamiltonian 
(|). If we define 

^-^(1 -i)- (10) 

the restrictions on the Hamiltonian (^|) are equivalent to the condition that the rotated matrix 
iUoHBdcUQ 1 is real antisymmetric, i.e. it belongs to the Lie algebra so(M), where M is the 
dimension of the Hilbert space (the same rotation of the Hamiltonian was used in || to identify 
the D universality class). 
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Figure 1: (a) Half-quantum vortex. Arrows denote the direction of vector d. (b) Spin vortex. The 
angle 9 is shown 

The last step in our argument is to note that, under the vortex boundary conditions, the 
dimension of the Hamiltonian (^|) is odd, not even (this may be difficult to visualize from the 
particle- hole representation ([)]), but easier from the rotated Hamiltonian UqH^qUq 1 ). Thus for 
a single- quantum vortex, we identify the space of the Hamiltonians as so(2N + 1) (class B in 
Cartan's notation, see Table 1). 



2 Half- Quantum Vortex 

Now let us consider a slightly different order parameter: suppose that the condensate wave function 
is described by the same expression (||), but now the vector d is able to rotate (either in plane or 
in all three dimensions) . The order parameter does not change under simultaneous change of sign 
of the vector d and phase shift by it of the phase ip: (ip, d) ~ {<p + ir, — d). This makes possible 
a half-quantum vortex with magnetic flux <I>o/2 On going around such a vortex, both the 

vector d and the phase ip rotate by ir (Fig. |]a) (see also discussion of such vortices in (l?], |l8|). 

Consider first the case when the vector d rotates in a plane (without loss of generality, let d 
rotate in the x-y plane). Then the condensate wave function (|J) is the sum of two decoupled 
components containing and [J.J,) pairings: 

(k x +iky) (11) 

[here r and 9 are the polar coordinates] . In a half-quantum vortex, one of the pairing components 
A^P and A? has an odd winding number of the phase around the vortex, while the other has an 
even winding number of the phase. In the simplest vortex of this type, the even winding number is 
zero, the corresponding component of the order parameter is nearly uniform and does not produce 
states near Fermi energy. The component with an odd winding number (without loss of generality, 
let it be A+\ and the winding number be one) has vortex boundary conditions. The mean-field 
Hamiltonian decouples and the vortex part of the Hamiltonian (spin-up sector) takes the form [we 



V(r,6)= A«(r,0)|tT> + Af(r,0)|U) 
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have also assumed an axially symmetric order parameter]: 

if = J d 2 r *{ - e F ^j * T + e i9 A(r)*| (V x + iV y ) + h.c 



In this form the half-quantum vortex is equivalent to a single-quantum vortex in a superconductor 
of spinless fermions (such a superconductor must have odd pairing). The quasiparticles do not 
have a definite spin projection, but mix particles and holes in the same spin-up sector: 7^ = 
/[u(r)^t(r) + w(r)$f (r)]<i 2 r. The Bogoliubov-deGennes equations for u and v turn out to be 
identical to those for a single-quantum vortex (^), and the solution of Kopnin and Salomaa [Of 
is equally applicable to this half-quantum vortex resulting in the same spectrum ([|). However, 
because of the relation rf{E) — j(—E), the number of subgap states in the half-quantum vortex 
is one half of those in a single-quantum vortex: particles in negative-energy levels are identical to 
holes in positive-energy levels. The zero-energy state is a Majorana fermion [jl8| 

Now add a disorder term to the Hamiltonian of a clean vortex dl3). Consider the following 
perturbations: 

• potential scattering: ffy = (r)^ a , 

• Zeeman field: ff h = (h(r)er a/3 ) ^p, 

• vector potential: ff a = ^ a (a(r) * p) 

• spin-orbit scattering: H\j = #3* (Uso(r) * [p X a a p}) */3, 

• deformation of order parameter: 

H A = t (A<g (r) * Px + A$ (r) * p„) *J + h.c. 

We assume that the perturbation may be switched on adiabatically in such a way that subgap 
levels stay localized. One possible way to fulfill this condition is to require that the perturbation 
is introduced only in a finite region around the vortex core. 

The symmetry classification proceeds slightly differently depending on whether the perturbation 
preserves the decoupling of the Hamiltonian into spin-up and spin-down sectors. 

The Hamiltonian remains decoupled under the perturbations Hy, H a , ffh with h || z, Hjj 
with Uso -L z, and if a with diagonal A a/ g. With these perturbations, the "spin- up" sector of the 
Hamiltonian in the basis (it, v) preserves its form (j^j with the same symmetries as for a single-flux 
vortex. The symmetry class of the Hamiltonian is therefore identified as so(2N + 1) (class B). 

If other perturbations are present, so that the Hamiltonian no longer splits into "spin- up" 
and "spin-down" parts, the full matrix of the Hamiltonian needs to be considered. Without any 
spin structure, the Hamiltonian breaking time-reversal and spin symmetries belongs to the so(M) 
symmetry class, as shown by Altland and Zirnbauer ||. Their argument also works in our case, 
with the reservation that the number of levels M is odd (with each level counted twice: once as 
a particle and the other time as a hole). It follows from the usual continuity consideration (the 
Majorana fermion at zero energy survives any local perturbation). Therefore the symmetry class 
is again so(2N + 1) (class B). [In the examples considered by Altland and Zirnbauer, the number 
of levels was even, which lead to the symmetry class D or so(2N)]. 

It is important however, that to reach the RMT limit in the latter case, such a system must 
have sufficiently strong disorder to bring the quasiparticles from the "spin-down" sector (gapped 
in the clean vortex) to the Fermi energy and to mix them strongly with the quasiparticles in the 
"spin-up" sector. In the present paper, we do not discuss conditions for reaching this limit. 
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3 Spin Vortex 



In this section we consider a hypothetical phase of a triplet superconductor analogous to the phase 
B of 3 He. Namely, the condensate wave function is assumed to be 



d*(k) ITT) + III) +id y (k)[ |TT) -III) 



d(k) = R 9 (k) , (13) 



and Re is a rotation in the plane by angle 9 (like in the previous sections, we consider a two- 
dimensional problem, and k is a two-dimensional vector). In this phase, the vector d rotates a full 
turn as the vector k goes around the circular Fermi surface. The parameter 9 denotes the angle 
between vectors k and d. We neglect spin-orbit interactions and assume that all values of 9 are 
allowed. 

It is important that this phase preserves the time-reversal symmetry. Moreover, a vortex 
structure is possible without time-reversal symmetry breaking. It is realized by rotating the angle 
9 a full turn of 2ir when going around the vortex center and keeping the phase <p constant (Fig. 0b). 
This vortex preserves the time-reversal symmetry and is a topological defect in the spin structure 
of the triplet pairing, thus we shall call it "spin vortex" . 

The Hamiltonian of the clean axially symmetric spin vortex splits into "spin-up" and "spin- 
down" sectors, similarly to the case of the half-quantum vortex considered in the previous section. 
In the spin vortex, both the "spin-up" and the "spin-down" components have vortex structure. The 



"spin- up" Hamiltonian is identical to the Hamiltonian H of (12), while the "spin-down" component 
is its complex conjugate H* (with spin reversed) — or vice versa. 

In a clean spin vortex, the spectrum of each of the two vortices (in the "spin-up" and the 
"spin down" sectors) are identical to that in the half-quantum vortex of the previous section, and 
therefore the total spectrum coincides with that of the single-flux vortex of Section I (with the two 
Majorana fermions combining into a single zero-energy level). In contrast to the single-quantum 
vortex, it is now not the spin symmetry which is responsible for the double degeneracy of the 
levels with non-zero energy, but the time-reversal symmetry (the degeneracy is due to Kramers' 
theorem). 

Now suppose that the Hamiltonian is perturbed by a disorder term preserving the time-reversal 
symmetry. The allowed perturbations include potential scattering Hy, spin-orbit scattering Hu, 
as well as deformation of the order parameter in a time-reversal invariant way. Then the levels 
with non-zero energy stay doubly degenerate, and the zero-energy level cannot shift (because it 
is non-degenerate). The symmetry class of the perturbed Hamiltonian may be identified from 
the argument of Altland and Zirnbauer as Dill (time-reversal symmetry present, spin symmetry 
broken) ||. In the p-wave vortex, the number of subgap levels is odd, and therefore the symmetry 
class is Dill-odd, in contrast to class Dill-even studied by Altland and Zirnbauer [similarly to 
distinguishing between classes B and D, we find it useful to distinguish between classes Dill-even 
and Dill-odd: they have different parameters a, as a consequence of the zero-energy level in odd 
dimensions! . 



4 Level Mixing by Disorder 

In this section we discuss whether the RMT limit may be realized in disordered vortices. We 
suppose that the superconductor contains a finite concentration ni mp of spinless impurities with 
strong electron scattering (Born parameter is of order one). Then the effect of impurities is 
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dual: they directly contribute a potential term Hy in the Hamiltonian and they also suppress 
the superconducting order parameter thus inducing an inhomogeneous pairing perturbation Ha- 
We shall see that the latter effect has a much stronger influence on shifting the energy levels of the 
subgap states. 

Consider first the potential term Hy with the impurity potential V(r) = £\ ViS(r — r^). The 
matrix element of the impurity potential between the two localized states is equal to 

{Hv)mn = ^2 v *P™n{ri) , /w(r) = «„(r)w„(r) - u^(r)u„(r)) . (14) 

i 

From explicitly solving Bogoliubov-deGennes equations for a clean p-wave vortex fll3| , one finds 
that u n and v n are proportional to each other, to the leading order in E/A (where E — uluq is 
the energy of the n-th level) . Therefore the matrix element p m n between any two low- lying states 
nearly vanishes (this was noticed by Volovik in |fl5f ) . A more accurate calculation shows that matrix 
elements of the charge density between states at energy E is of the order of £~ 2 (E/ A), which is 
by the factor (E/A) smaller than that in the s-wave vortex. As a consequence, even for strong 
impurities (with Vi ~ £f kp 2 ), an extremely high impurity concentration (rti mp ~ k 2 F ) is required 
to mix low-energy levels — a much higher disorder than that destroying p-wave superconductivity 
("■imp ~ &f£ -1 ) [ ^oj |2l[l . Therefore, for a uniformly disordered p-wave superconductor, the Hy 
term is insufficient for the RMT regime in the vortex core. 

The suppression of the order parameter H A turns out to be a more important effect. In 
unconventional superconductors, spinless impurities suppress pairing in the region of size £ which 
is of the same order as the size of the vortex core. For a homogeneous order parameter (without 
vortex), the suppression of the order parameter <5A(r) decreases as (in our two-dimensional 
problem) at distances up to £, and is negligible at distances much larger than £ [^2[ and the 
integral suppression may be estimated as J 5A d 2 r ~ A^kZ . Also, in a chiral superconductor, the 
component of the order parameter of opposite chirality is admixed in the vicinity of impurity p3 . 
The two effects gives comparable contributions to the matrix elements between subgap states. An 
accurate calculation of the impurity influence on the states in the vortex core requires a full self- 
consistent solution for the order parameter taking into account both the vortex and the impurities. 
This goes beyond the scope of this paper, and we only estimate the order of magnitude of the 
corresponding matrix elements. The deformation of the order parameter Ha has the matrix 
elements between subgap states m and n: 

(H A ) mn = Jd 2 r ( u* m (r) [<5A^(r) * (*V a )] u n (r) 

+v* m (r) [sA^*(r) * (<V a )] v n (r)j . (15) 

Unlike the case of the matrix elements of Hy, there is no cancellation in {H^,) mn , and we may 
estimate, for one impurity in the vortex core 



(H A ) mn ~ 1 J 5A(v) d 2 r ~ A(fc i ,£)- 1 



^0 • (16) 



Thus, a single impurity introduces matrix elements of the order of interlevel spacing. These 
estimates suggest that the RMT regime may be achieved when the number of impurities in the 
vortex core is much greater than one (which corresponds to the moderately clean regime A -1 <C 



5 Summary and Discussion 



We have shown that symmetries of vortex Hamiltonians in p-wave superconductors correspond 
to one of the two random-matrix ensembles: B or Dill-odd. These two RMT ensembles are 
distinguished from the rest of ensembles by the presence of zero-energy modes (zero-energy modes 
are also present in chiral ensembles at p ^ q, see Table |I|). Zero-energy levels appear in p- 
wave vortices as a consequence of odd pairing symmetry combined with the vortex (topologically 
nontrivial) boundary conditions. We have checked the symmetries of the Hamiltonian explicitly 
taking three vortices as examples: a single-quantum vortex with spin symmetry, a half-quantum 
vortex, and a flux-less (spin) vortex without time-reversal symmetry breaking. Based on these 
examples, we may conjecture that classes B and Dill-odd appear in any vortex-like structure with 
odd pairing, whenever a zero-energy mode is present. In cases with time reversal symmetry, the 
vortex belongs to the class Dill-odd; when the time-reversal symmetry is broken — to the class 
B. Of course, this can also be checked explicitly in any particular case. 

The two classes B and Z?III-odd considered in this paper are odd-dimensional counterparts of 
the classes D and Dill-even, respectively (see Table 1), with the even-dimensional classes realized 
in conventional (singlet-pairing) normal-superconducting structures || . Far from zero energy, the 
statistics of energy levels is not affected by the zero mode. In the immediate vicinity of zero energy, 
repulsion from the zero mode in classes B and Dill-odd increases the exponent a (see Table 1), 
suppressing the average density of states (p(u))) as (p(w)) <x w 2 for class B, and {p(u>)) oc \uj\ 5 for 
class Z?III-odd. 

Because of the zero modes, the classes B and -Dill-odd stand somewhat separately from the 
rest of RMT ensembles. In particular, they do not appear in the table of correspondence between 
symmetries of the Hamiltonian and those of the transfer matrix || |(| (for RMT ensembles without 
zero-energy levels, this table establishes a one-to-one correspondence). Neither these two ensembles 
are known to be derived from a supersymmetric sigma- model as other classes are Q . Understanding 
the latter fact remains an interesting problem, as well as a possible microscopic derivation of the 
random-matrix theory for p-wave vortices with the supersymmetric method, analogously to that 
for the s-wave vortex . 

Several more comments can be made on comparing our results to the level statistics in s-wave 
vortices. The latter is known to belong to the symmetry class C [^5[ [26) having spin-rotational 
symmetry and broken time-reversal symmetry ||. In contrast with the SU(2) spin symmetry 
in conventional superconductors, the single-quantum vortex considered in our paper has the 0(2) 
spin symmetry instead, which leads to a somewhat different symmetry classification. The resulting 
symmetry is so(iV), which leads to classes D (in even dimensions, no vortex) or B (odd dimensions, 
vortex with zero mode). 

The joint probability distributions of the energy levels (||) for classes B (p-wave vortex) and C 
(s-wave vortex) coincide (with a = /3 = 2), and only the zero-energy level distinguishes the two 
ensembles. In particular, the average density of states in the class B random-matrix ensemble is 

/ f \\ 1 sin(27ro;/tj ) 

(p(u))) = +6(uj). (17) 

LUo 2,71 LO 

(for class C it is the same but without the 6(u>) term). 

In s-wave vortices, Koulakov and Larkin found that in a wide range of impurity concentrations, 
the level mixing does not lead to a class C, but has more symmetries |^7j . These symmetries arise 
from the fact that the impurities are local scatterers in a chiral quasi-one-dimensional two-channel 
system, resulting in an ensemble of 2 x 2 random matrices, instead of N x N matrices with large 
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N. Such effect seems improbabfe in p-wave vortices where the main contribution to ievel mixing 
comes from the distortion of the order parameter around impurities which is extended in space 
to distances of the order of the core size. In this paper we performed only an order-of-magnitude 
estimate suggesting that the RMT limit may be reached in moderately clean vortices (with the 
number of impurities per vortex core much greater than one). A more accurate self-consistent 
treatment of impurities in the vortex core is required for a quantitative study of the effects of 
disorder. 

The author wishes to thank M. V. Feigel'man for suggesting this problem, discussions, and com- 
ments on the manuscript. At different stages of work, the author benefited from discussions with 
G. Blatter, V. Geshkenbein, D. Gorokhov, R. Heeb, L. Ioffe, C. Mudry, M. Skvortsov, G. E. Volovik, 
and M. Zhitomirsky. The author thanks Swiss National Foundation for financial support. 
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